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Using a lGW-lps pump laser pulse in high gain parametric down-conversion allows us to detect 
sub-shot-noise spatial quantum correlation with up to one hundred photoelectrons per mode, by 
means of a high efficiency CCD. The statistics is performed in single-shot over independent spatial 
replica of the system. The paper highlights the evidence of quantum correlation between symmetrical 
signal and idler spatial areas in the far field, in the high gain regime. In accordance with the 
predictions of numerical calculations the observed transition from the quantum to the classical 
regime is interpreted as a consequence of the narrowing of the down-converted beams in the very 
high gain regime. 

PACS numbers: 



Spatial quantum optical fluctuations are studied be- 
cause of new potential applications of quantum optical 
procedures in parallel processing and multi-channel op- 
eration. Examples are quantum holography , the quan- 
tum teleportation of optical images [2| , and the measure- 
ments of small displacements beyond the Rayleigh limit 
There is now a large literature on spatial effects in 
the spontaneous regime of parametric down-conversion 
(PDC) where photons are created one pair at a time Q. 
The process of PDC is in fact particularly suitable for the 
study of spatial correlations because of its large emission 
bandwidth in the spatial frequency domain [5(. Never- 
theless to date most spatial correlation measurements in 
PDC have been performed in single photon counts regime 
0, without evidencing any relevant quantum effects. 
The quantum twin beam character of the PDC emission 
has been evidenced in || by using low-to-medium pump- 
power lasers (< 1 MW) and relying on statistical ensem- 
bles from different temporal replicas of the system. With 
increasing gain a transition from the quantum to the clas- 
sical regime has been observed 0- However, recent the- 
oretical investigations predict multi-mode spatial quan- 
tum correlations (sub-shot-noise photon-number correla- 
tion between symmetrical portions of the signal and idler 
angular distributions) also for arbitrarily high gains, pro- 
vided that the detection area exceeds the typical size of 
the mode (coherence area) 0, ^) . 

Here we report on the first quantum spatial mea- 
surements of PDC radiation performed by using a low- 
repetition rate (2 Hz) pulsed high-power laser (1GW- 
lps). This enables us to tune the PDC to the high-gain 
regime while keeping a large pump beam size (~ 1mm). 
The huge number of transverse modes (roughly given by 
the ratio between: (i) the area of the near-field gain pro- 
file and (ii) the inverse of the angular bandwidth of the 
PDC process) allows us to identify regions of the sig- 
nal and idler beams where symmetrical signal-idler pixel 
pairs correspond to independent spatial replica of the 



quantum system. We concentrate on a portion of the 
parametric fluorescence emitted close to the collinear di- 
rection and within a narrow frequency bandwidth around 
degeneracy The generated pairs of signal and idler 
phase-conjugate modes propagate at symmetrical angles 
with respect to the pump direction in order to fulfil the 
phase-matching constraints, and each pair of symmetri- 
cal spots characterizing the far field represents a spatial 
replica. Thanks to the very large number of these, the 
statistical ensemble averaging necessary for the quantum 
measurement can be solely done over the spatial replicas 
for each, single, pump-laser pulse. A characterization of 
the system over several shots was only made in our case 
in order to verify that the selected spatial replicas are 
indeed statistically identical, as required for the suitable 
definition of the ensemble. The single-shot measurements 
reveal sub-shot-noise spatial correlations for a PDC gain 
G (intensity amplification factor) leading to the detec- 
tion of up to ~ 100 photoelectrons per mode. Finally, by 
numerically solving the three-waves coupled equations in 
the framework of a 3D+1 quantum model, we are able to 
attribute the observed transition from quantum to classi- 
cal regime to the near-field gain narrowing that occurs in 
presence of a bell-shaped pump beam, at very-high gain. 

The experimental setup is sketched in Fig. ^ The 
third harmonic (352 nm) of a lps, chirped-pulse ampli- 
fied Nd:glass laser (TWINKLE, Light Conversion Ltd.) 
is used to pump a type II 5x7x4mm 3 /3-barium borate 
(BBO) non-linear crystal, operated in the regime of para- 
metric amplification of the vacuum-state fluctuations. 
The input and output facets of the crystal are anti- 
reflection coated at 352 nm and 704 nm, respectively. 
The pump beam is spatially filtered and collimated to a 
beam waist of approximately 1 mm (FWHM) at the crys- 
tal input facet. The energy of the 352nm pump pulse can 
be continuously tuned in the range 0.1-0.4 mJ by means 
of suitable attenuating filters and by changing the en- 
ergy of the 1055nm pump laser pulse, allowing to have 



a gain in the range 10 < G < 10 3 . The parametric 
fluorescence at the horizontally polarized signal and ver- 
tically polarized idler modes is emitted over two cones, 
whose apertures depend on the specific wavelengths (see, 
e-g-, Q)- The BBO crystal (6 = 49.05°, cf> = 0) is 
oriented in order to generate signal and idler radiation 
cones tangent in the collinear direction at the degener- 
ate wavelength lu s = uji = uj p /2 (s, i and p referring to 
signal, idler and pump respectively). The fluorescence 
around the collinear direction is selected by means of a 
5mm x 8mm aperture, placed 15 cm from the output 
facet of the BBO. The aperture turned out to be nec- 
essary in order to prevent beam clipping by the PBS, 
otherwise giving rise to substantial scattered radiation. 
The selected portion of the beam is transmitted through 
a polarizing beam splitter (PBS), which separates the 
signal and idler beams. The latter are finally sent onto 
two separate regions of a deep-depletion back illuminated 
charged coupled device (CCD) camera 0] (Roper Scien- 
tific NTE/CCD-400EHRBG1, 16 bits dynamical range, 
quantum efficiency r)K, 89% at 704 nm at -40°C, dark cur- 
rent and read out noise < le/pixel/s and < 3e/pixel/s 
respectively), placed in the common focal plane of the 
two lenses (/=10 cm) used to image the signal and idler 
far fields. The detection array has 1340x400 pixels, with 
a pixel size of 20/imx20^m. Prior to the experiment the 
CCD was calibrated with a coherent source allowing the 
retrieval of spatial shot-noise statistics in its full dynamic 
range [l6| . In our setup the correlation measurements are 
performed without using any narrow-band interferential 
filters (IFs), in contrast to the case of photon-counting 
experiments (coincidence measurements), since IFs un- 
avoidably introduce relevant transmission losses reducing 
the visibility of sub-shot noise correlations. The pump- 
frequency contribution is removed by using normal in- 
cidence (M 5 ) and at 45 deg (M4) high-reflectivity (HR) 
mirrors coated for 352 nm placed before and after the 
PBS, respectively, and a low-band pass color filter (90% 
transmission around 704 nm) placed in front of the CCD. 
A further HR@352nm mirror (M4) is placed in one of 
the two arms at a suitable angle in order to balance the 
unequal transmission of the PBS in the two arms. All 
the optical components (except the color filter) have an- 
tireflection coatings at 704 nm. The estimated quan- 
tum efficiency of each detection line, which accounts for 
both the transmission losses and the detector efficiency, 
is ritot ^75%. 

Fig. shows a typical far field image recorded in a 
single shot, where a fairly broadband radiation (i.e. the 
one transmitted by the rectangular aperture) is acquired 
in the signal (left) and idler (right) branches. The se- 
lection of the desired temporal and angular bandwidth 
around degeneracy is made by temporarily inserting in 
front of the CCD a lOnm wide IF around 704 nm, allow- 
ing us to locate the collinear degeneracy point (see Fig. 
|2h>) ■ The data analysis is limited within two rectangu- 
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FIG. 1: Scheme of the experimental set-up used for the spatial 
correlation measurements (see text). 



lar boxes (black frames in Fig[2^) corresponding to an 
angular bandwidth of 20 mrad x 8 mrad and to a tempo- 
ral bandwidth smaller than 10 nm. The selected regions 
contain 4000 pixels each. In this work we investigate 
pixel pair correlation, and since the size of the CCD pixel 
approximately corresponds to the physical size of the 
replica, the ensemble is large enough to perform the de- 
sired statistics. We have observed that much larger boxes 
worsen the level of the correlation; this is attributed to 
residual scattering owing to diffraction from the borders 
of the aperture but also to the contribution of other fre- 
quency components (far from degeneracy) characterized 
by a lower angular symmetry between signal and idler 
cones. A zoom of the selected areas is presented in Fig(2t ? 
where the rather spectacular symmetry of the intensity 
distribution in the signal and idler branches shows the 
twin-beam character of the phase-conjugate modes. 

The aim of the experiment is to quantify the sym- 
metrical pixel pair correlation. This is done by mea- 
suring the variance <J 2 _ i of the PDC photoelectron- 
number difference n s — rii of the signal/idler pixel pair 
versus the mean total number of down-converted pho- 
toelectrons (pe) of the pixel pair. This variance is 



{(n s - rii) 2 ) - {n s 



i) 2 where the averages are 



spatial averages performed over all the symmetrical pixel 
pairs contained in the chosen regions. Each single shot 
of the laser provides a different ensemble, characterized 
by its pixel pair average pe number (n s + in turn 
related to the parametric gain. In the experiment, en- 
sembles corresponding to different gains are obtained by 
varying the pump-pulse energy. We note that the read- 
out noise of the detector, its dark current, and some un- 
avoidable light scattered from the pump, signal and idler 
fields contribute with a non-negligible background noise 
to the process. This is taken into account by applying a 
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FIG. 2: (a) Single-shot far field image recorded by the CCD 
for a pump beam waist wo ~ 1mm and pump energy e p ~ 
0.3mJ. The spatial areas for statistics are delimited by the 
white boxes selected within the degenerate signal and idler 
modes, spatially localized from the single shot image recorded 
with the lOnm-broad IF (b). (c) Zoom of two symmetrical 
areas of the signal and idler far fields. 



standard correction procedure (see for example [18j), by 
subtracting the background fluctuations a\ from the ef- 
fectively measured variance cr? s+b \, i+f) , of the total inten- 
sity difference (signal+background)- (idler +background) 
obtaining a 2 s _ i = o? a+b \_/ i+b \ — ^ a b- This background 
noise, having a standard deviation of 7 counts (±0.1 from 
shot to shot, estimated by repeating the measurement 
with the same pump-pulse energy) is measured in pres- 
ence of pulse illumination over an area of the same size 
of the acquisition area and suitably displaced from the 
directly illuminated region. The validity of the data cor- 
rection procedure is tested by sending in the setup (with 
no crystal) through the PBS a coherent circularly polar- 
ized pulsed beam (@704nm), and verifying for different 
laser energies that the intensity difference fluctuations 
from the two coherent portions of beams recorded on the 
CCD lie at the shot noise level. 

Fig-Glshows the experimental results where each point 
is associated with a different laser shot. The data are 
normalized to the shot noise level (SNL), and their statis- 
tical spread accounts for the background correction. Al- 



FIG. 3: Intensity difference variance o 1 s _ i normalized to the 
SNL (n s + m)- Each point (white circle) corresponds to a 
single shot measurement where the spatial ensemble statis- 
tics has been performed over a 100 x 40 pixels region. The 
triangles (each one obtained by averaging the experimental 
points corresponding to a certain gain) and their linear fit il- 
lustrate the trend of the data in the region between (n s + rii) 
=8 and 20. Inset: Typical correlation degree profile in the 
regime where (n s + m) ~ 8 (see text). 

though the noise on the individual signal and idler beams 
is found to be very high and much greater than the stan- 
dard quantum limit (={n s ) and (rij) respectively), we 
observe an evident sub-shot noise pixel pair correlation 
up to gains characterized by (n s + rii}~ 15 — 18. Since 
in that regime the observed transverse size of the co- 
herence areas (i.e. of the modes) is about 2-4 pixels, 
this approximately corresponds to 100 pe per mode. We 
can have an idea of the transverse size of the mode by 
looking at the standard two-dimensional cross-correlation 
degree 7 = ((n s rii) — (n s )(ni))/ \J cr^crf, between all the 
angularly symmetrical signal and idler pixels contained 
within the black boxes (see Fig. |2Ji) . This can be plotted 
for instance as a function of the horizontal and vertical 
shifts of the recorded image on the CCD, keeping fixed 
the position of the boxes. In general | 7 |< 1 with 7 = 1 
for perfect correlation. A transverse section of the corre- 
lation function obtained from a single-shot image char- 
acterized by {n s + 7ij}« 8 is plotted in the inset of Fig. 
121 as a function of the horizontal shift x (in pixels units). 
As expected, virtually perfect correlation (in our case the 
peak value is ~ 0.99) is obtained for perfect determina- 
tion (i.e. within one pixel) of the center of symmetry 
between the signal and the idler regions. 

In order to interpret the observed transition from quan- 
tum to classical regime we present in Fig. 0] the results of 
the numerical calculations. The full quantum model ac- 
counts for the two transverse and the temporal degrees 
of freedom with propagation along the crystal, for the 
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FIG. 4: Numerical calculation of u 2 a _ i (normalized to SNL) 
between symmetrical portions of signal and idler plotted as a 
function of the detection area represented by N x N binned 
pixels. Different curves correspond to different values of the 
gain characterized by the mean number of down-converted pe 
per pixel pair (n s + n,i). 



angular and chromatic material dispersion up to the sec- 
ond order, and for the finite spatial and temporal widths 
of the Gaussian pump pulse. The crystal and input- 
pulse parameters are those relative to the experiment. 
The figure presents <rf_j, normalized to the SNL, vs the 
size of the detection area for different gains. Each point 
is the result of a statistics performed over one single 
laser-shot. The case N=l corresponds to the experiment. 
The simulations (data not shown) outline that, in spite 
of the fixed pump-beam diameter, the signal and idler 
beam diameters at the crystal output strongly depend 
on the gain and decrease when the latter increases. This 
can be easily interpreted when considering that the sig- 
nal and idler beam size maps not the pump-beam pro- 
file but the actual parametric amplification gain profile 
G(r) ~ cosh 2 [erA(r)L] (L being the crystal length, A 
the pump field amplitude and a a parameter proportional 
to the setting characteristics), as long as filtering due to 
the limited spatial bandwidth does not take place [2fij |. 
On narrowing the size of the PDC beams, the coherence 
areas in the far field (i.e. the modes) increase their size, 
as straightforward consequence of the convolution theo- 
rem in Fourier analysis |l3| . Since revealing quantum cor- 
relations requires detection areas lar ger (or comparable) 
to the mode size (as also discussed in )l(|). it is necessary 
when increasing the gain to have larger detectors in or- 
der to obtain below-shot-noise variance as shown in Fig. 
QJ Note that Fig. 0] evidences the transition from quan- 
tum to classical regime in case of single-pixel detection 
(N=l) for a gain that is higher than in the experiment. 
Instead, in the experiment, excess noise is observed for 
(n s + rii) >20, which we attribute first to the effect of 
residual scattered light whose contribution grows linearly 



with the radiation fluence and is thus expected to over- 
come the shot noise at large pumping, and second to the 
uncertainty in the determination of the symmetry center 
of the signal and idler image portions. In fact simulations 
have shown that an uncertainty as small as a few microns 
(i.e. a fraction of the pixel size, unavoidable experimen- 
tally), prevents to observe sub-shot-noise correlation as 
soon as (n s + rn) exceeds some tens of pe, while still pre- 
serving sub-shot-noise correlation for smaller gain values. 
Finally, the maximum level of noise reduction observed 
experimentally agrees with the theoretical limit (dotted 
line in Fig-El determined by the total losses of the system 
(~ 1 — r\tot ,10]), in accordance with the result of Fig. 0] 
In conclusion, we have shown that twin beams of light 
generated in parametric down-conversion exhibit sub- 
shot noise spatial correlation by measuring an evident 
quantum noise reduction on the signal/idler intensity dif- 
ference. A transition to above shot-noise correlation is 
observed as the gain increases. This quantum-to-classical 
transition, in agreement with numerical simulations, is 
explained as a narrowing of the signal/idler beams with 
increased gain. This leads in turn to a larger mode size 
and therefore also to the need of larger pixels to observe 
below shot-noise correlation [Jjj . This will be the aim of 
a future work. To our knowledge, this is the first exper- 
imental investigation of quantum spatial correlations in 
the high gain regime, where the huge number of trans- 
verse spatial modes is detected in single shot by means 
of a high-quantum-efficiency CCD. 
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